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Abstract 

Assuming geometric Brownian motion as unaffected price process S^, Gather al & 
Schied (2011) derived a strategy for optimal order execution that reacts in a sensible 
manner on market changes but can still be computed in closed form. Here we will in- 
vestigate the robustness of this strategy with respect to misspecification of the law of 
S^. We prove the surprising result that the strategy remains optimal whenever is a 
square-integrable martingale. We then analyze the optimization criterion of Gatheral &; 
Schied (2011) in the case in which is any square-integrable semimartingale and we 
give a closed-form solution to this problem. As a corollary, we find an explicit solution to 
the problem of minimizing the expected liquidation costs when the unaffected price pro- 
cess is a square-integrable semimartingale. The solutions to our problems are found by 
stochastically solving a finite-fuel control problem without assumptions of Markovianity. 

Key words: market impact, optimal order execution, Almgren-Chriss model, robustness, 
model uncertainty 



1 Introduction 

This paper can be read from two complementary perspectives. 

From the first perspective, it is a paper on the optimal execution of large orders, which 
is a problem that was first discussed by Bertsimas &: Lo (1998) and Almgren & Chriss 
(1999, 2000). The construction of optimal order execution strategies has considerable 
practical significance. As always, strategies to be applied in practice should have reason- 
able quantitative and qualitative properties, should be easy to implement, and, ideally, 
robust with respect to model misspecification. The first properties are satisfied by the 
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strategy that was derived in Gatheral & Schied (2011). Assuming geometric Brownian 
motion as unaffected price process S^, this strategy was obtained as the minimizer of a 
cost functional, which can be regarded as the time-averaged risk of the remaining posi- 
tion. In this paper, we investigate how the optimality of this strategy (x^) is affected 
by changes of the distribution of S^. To this end, we will give a closed- form solution of 
the strategy that minimizes the cost criterion from Gatheral & Schied (2011) if is a 
general square-integrable semimartingale. Surprisingly, it will turn out that the optimal 
strategy coincides with (x^) whenever is a square-integrable martingale. In this sense, 
(x^) is very robust with respect to misspecification of the dynamics of the unaffected price 
process and, as a consequence, satisfies most requirements one would have on a reasonable 
order execution strategy. 

This brings us to the second perspective, from which this paper can be viewed as a case 
study in robustness with respect to model uncertainty. The ubiquitous existence of model 
uncertainty was first emphasized by Knight (1921), but only few systematic approaches 
to this phenomenon exist to date; we refer to Cont (2006) and, for an overview over some 
material, to Follmer & Schied (2011) and the references therein. The present paper adds a 
further particular to the literature on model uncertainty, namely the study of robustness 
with respect to model misspecification in a stochastic control problem arising in order 
execution. This feature of robustness is closely related to the remarkable fact that here 
it is possible to solve explicitly a stochastic control problem with fuel constraint without 
assumptions of Markovianity and without using partial differential equations. 

As a corollary of our main results, we are able to give an explicit solution to the 
problem of minimizing the expected liquidation costs in the Almgren-Chriss framework 
when the unaffected price process is a square-integrable semimartingale. We find that 
optimal strategies always exist and that the drift enters the corresponding formula in 
integrated form. One can thus expect that possible misspecifications of the drift may 
average out. This relatively stable behavior is in stark contrast to the direct dependence 
of optimal strategies on the derivative of the drift in models with transient price impact 
as found in Lorenz & Schied (2012). 

In the subsequent Section 2 we will explain in some detail the background for our study 
and we will provide a precise formulation of the problem we are looking at. In Section 
3 we will state our main results. We will start by formulating the results pertaining to 
martingale dynamics of the unaffected price process. These results are, however, just 
corollaries of our main result. Theorem 2, which provides the closed- form solution of the 
optimal strategy general semimartingale dynamics of S^. As a further corollary, we find 
an explicit solution to the problem of minimizing the expected liquidation costs when 
the unaffected price process is a square-integrable semimartingale. Our formula for the 
optimal strategy in Theorem 2 is obtained by first guessing the minimal value of the 
optimization problem and then applying a stochastic verification argument to confirm the 
guess. In Section 4, we have included heuristic arguments, based on partial differential 
equations, which show how this guess can be found. All proofs are provided in Section 5. 

2 Background and problem formulation 

In the continuous-time version of the market impact model of Almgren & Chriss (1999, 
2000) it is assumed that the number of shares in the portfolio of a trader is described by 
an absolutely continuous trajectory t i— >• xt- Given this trading trajectory, the price at 
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which transactions occur is 



= S'^t +r]Xt + -i{xt-xo), (1) 

where > and 7 > are constants and S'f is the unaffected stock price process. The 
term r]Xt corresponds to the temporary or instantaneous impact of trading Xt dt shares 
at time t and affects only this current order. The term ^{xt — xq) corresponds to the 
permanent price impact that has been accumulated by all transactions until time t. The 
unaffected price process is usually assumed to be a martingale. There are good reasons, 
however, why it can make sense to relax the martingale assumption. For instance, there 
may be other large traders active in the market and their trading activities create a drift 
on top of random market fluctuations; see e.g., Schoneborn &; Schied (2009). 

Let us now consider an order execution strategy in which an initial long or short 
position of X shares is liquidated by time T. The asset position of the trader, {xt)Q<t<T^ 
thus satisfies the boundary condition xq = X and xt = 0. In such a strategy, —xt dt 
shares are sold at price 5f at each time t. Thus, the costs arising from the strategy 



{xt)o<t<T are 

C{x) := [ 
Jo 



T 

Sfxt dt 



T 



--XSo- [ xtdS^ + rj [ ±fdt + lx^, 
Jo Jo ^ 

where we have used integration by parts. 

The optimal order execution problem consists in maximizing a certain objective func- 
tion, which may involve revenues and additional risk terms, over a suitable class of admissi- 
ble trading strategies {xt)o<t<T with side conditions xq = X and xt = 0. The easiest case 
corresponds to minimizing the expected costs when is a martingale. In this case, the 
expectation of the stochastic integral xt dS^ vanishes for suitably bounded strategies, 
and we obtain 

T 



E[C(x)] = -xo5o + ^xg + r/Ef / xldt\ (3) 
^ ^ Jo ^ 





In this setting, minimization of the expected costs was first considered in Bertsimas & Lo 
(1998) in a discrete-time framework. A simple application of Jensen's inequality shows 
that the unique strategy that minimizes the expected costs (3) is characterized by having 
the constant trading rate 

■VWAP _ _xo 

regardless of the particular dynamics of the martingale S^. When, as is usually assumed in 
practice, time is parameterized in volume time, such a constant trading rate corresponds 
to a VWAP strategy, where VWAP stands for volume-weighted average price. 

First problem: Minimize the expected costs E[C(x)] when is not a martingale but 
a general square-integrable semimartingale. 

Surprisingly, the preceding problem can be solved explicitly in full generality. Our cor- 
responding result. Corollary 2, will be derived as a special case of a more general result, 
Theorem 2. To motivate its statement, we need to look into cost-risk criteria that go 
beyond the expected costs of an order execution strategy. 
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Almgren & Chriss (1999, 2000) were the first to point out that executing orders late 
in the trading interval [0, T] incurs volatihty risk. They therefore suggested to minimize 
a mean-variance functional of the form 

E[C(x)] + |var(C(x)), (4) 

where a is a risk-aversion parameter. While mean-variance optimization may be ap- 
pealing to practitioners due to its common use in finance, this approach has two major 
disadvantages when applied in order execution. First, it is not easy to find mean- variance 
minimizing strategies unless one restricts strategies to be deterministic and assumes that 
the unaffected price process is a Bachelier model, 

= Sq + aWt with a ^ and W a Brownian motion. (5) 

In this latter case, calculus of variations easily yields 



as the unique deterministic mean-variance optimal strategy. The second disadvantage 
stems from the fact that the mean- variance functional (4) is not time- consistent, since it 
involves a squared expectation operator. As a consequence, an optimal adaptive strategy 
computed at time t = loses its optimality at any later time, even if market conditions 
remain unchanged. Another consequence of time inconsistency is that techniques from 
stochastic optimal control cannot be applied directly, which greatly complicates the com- 
putation of mean-variance minimizing strategies. We refer to Lorenz & Almgren (2010) 
and Forsyth (2011). 

The time consistency of the optimization problem can be retained for the maximization 
of expected utility, 

E[U{-C{x))], 

where ?7 : R — )■ M is a concave, increasing utility function. The strategy maximizing 
the expected utility in the class of all adaptive strategies can be characterized by means 
of a nonlinear Hamilton-Jacobi-Bellman (HJB) partial differential equation (PDE) with 
singular initial condition; see Schied &: Schoneborn (2009) or Schoneborn (2011). This 
equation usually cannot be solved in explicit form unless U{x) = —e~^^ and assumption 
(5) holds, in which case we recover (6) as optimal strategy; see Schied et al. (2010). In all 
other cases, numerical techniques for solving nonlinear PDEs with singular initial condi- 
tion will be necessary. Moreover, optimal strategies may have counterintuitive behavior 
in reaction to certain parameter changes; see (Schied & Schoneborn 2009, pp. 190-191). 

We also refer to Forsyth et al. (2009) for another risk criterion that also leads to a 
singular HJB equation, which is similar to the one found in the maximization of expected 
utility. 

To summarize, all of the optimization criteria we have discussed so far have at least 
one of the following four disadvantages: 

• they yield only deterministic strategies that do not react on the movement of asset 
prices; 

• they are time-inconsistent; 
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• their computation requires complex numerics for solving a nonlinear PDE with sin- 
gular initial condition; 

• or they admit counterintuitive behavior in reaction to certain parameter changes. 

These properties are all not desirable from a practical point of view. Gatheral & Schied 
(2011) therefore proposed another optimization criterion, which leads to a strategy that 
is sensitive to changes in the asset price, that can be easily computed in closed form, 
and whose reaction to parameter changes is completely transparent. This optimization 
criterion is based on the common practice in risk management to assess the risk of a 
position of X > shares as a constant multiple of their current value. Thus, the risk 
of the asset position xt at time t is assessed as XxtSf for some constant A > 0. This 
constant A is typically derived from the Value at Risk of a unit asset position under the 
assumption of log-normal future returns. As argued in (Gatheral & Schied 2011, Remark 
2.2), one could obtain the same formula (but perhaps with a different value of A) if Value 
at Risk is replaced by a coherent risk measure or by any other positively homogeneous 
risk measure. The optimization criterion proposed in Gatheral & Schied (2011) consists in 
minimizing the following sum of the expected execution costs and the expectation of the 
time-averaged Value at Risk of the positions xt held during an order execution strategy 



E 



C{x)+ / XxtSfdt 







To simplify notations, we will henceforth consider the minimization of the cost functional 

rT 



E 



C{x)+ / \xt{S^ + jxt) dt 







(7) 



which can be easily transformed into the minimization of E[C(x) + Jq XxtSf dt]. 

In Gatheral & Schied (2011), the minimization of the functional (7) was considered 
in the case where the unaffected price process (5^) is a risk- neutral geometric Brownian 
motion, 

50 = 5oe"^*-3-'*. (8) 

It was stated in (Gatheral & Schied 2011, Theorem 3.2) that for 7 > the optimal 
admissible strategy I 
is given by 



minimizing the functional (7) within a suitable class of strategies 



sinh {u{T - t)) 



X 



sinh {uT) 2u Jq 1 + cosh {v{T - s)) 



ds 



X/t] and z/^ = X^/rj. The optimal strategy in the limiting case 7 



(9) 



is 



t 



T 



X 



5? ds 



(10) 



where A : 
given by 

AT 

~J0 

see (Gatheral &: Schied 2011, Theorem 3.1). 

The solutions (x^) and (x^) are clearly adaptive and react on changes of the asset price. 
More precisely, they are aggressive in the money in the sense that shares are sold faster 
when stock prices go up. Moreover, in the case 7 > we find for the choice A = aa'^ /2'y 
that 

^*_^MV Asinh(i/(r-t)) 

/o 



2u 



1 + cosh {u{T - s)) 



ds, 
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where is as in (6) . So (x| ) liquidates a given asset position fast than the deterministic 
mean- variance optimal strategy (xj^^). 

A remarkable property of the optimal strategies (x^) and (x^) is that they are indepen- 
dent of the volatility a of the unaffected price process . This feature already indicates 
a certain robustness of (xj') and (x^) with respect to model uncertainty. The second goal 
of this paper is to analyze the robustness of (xj') and (x^) in a systematic manner: 

Second problem: In setting up our optimization problem (7), we have assumed that 

follows a geometric Brownian motion. But suppose that in reality has different 
dynamics. How will this affect the optimality of our strategies (x^) and (x^)? 

In the next section, we will approach this problem in two steps. In the first step, we 
will assume that is a martingale. In this case, the surprising answer to our question 
will be that the strategies (x^) and (xj ) remain optimal whenever is a rightcontinuous 
and square- integrable martingale. In this sense, the optimal strategies {xD and (x^) are 
very robust; their optimality depends only on the martingale property and not on the 
particular distribution of . In the next step, we will also drop the martingale property 
and assume only that is a square-integrable semimartingale. By formally taking A = 0, 
we will then obtain the solution to our first problem as a special case. 



3 Main results 

Let us start by formally setting up the optimization problem. All stochastic processes shall 
be defined on a filtered probability space (1^,7-", (J-t),P) that satisfies the usual conditions 
and for which Tq is P-trivial, i.e., T[A] € {0,1} for all A ^ Tq. The unaffected price 
process is assumed to be a cadlag semimartingale that is square-integrable in the sense 
that 

E[( sup \S^^\f] <oo. (11) 

0<t<T 

By X{T, X) we denote the class of all admissible strategies for the problem of liqui- 
dating X > Q shares during the time interval [0, T] (analogous statements hold for the 
problem of buying a position of X > shares). This class consists of all adapted and 
absolutely continuous strategies that satisfy the side conditions xq = X and xt = and 
the integrability condition 



E 



T 



X? dt 







< oo. (12) 



This condition is clearly necessary for (7) to make sense and to be finite. In fact, it is also 
sufficient: 



Proposition 1. Under assumption (11), the cost functional (7) is well-defined and finite 
for any x € ?(!{T,X). If, moreover, is a martingale, then the following identity holds: 



E 



T 



C(x) + / XxtiS° + 7x0 dt^ = - XSo + 7?E 



[xt + XSfxt + ly^x^) dt 



where again X = X/rj and = X^/rj. 



6 



We can now state our first result. 

Theorem 1. Assume that is any rightcontinuous and square-integrable martingale 
satisfying (11). For v = Xj/r] > 0, the unique strategy minimizing the cost functional 



E 



T 



(xt + XStXt + u^xf) dt 



^ JO 



within X{T,X) is given by (9) and the value of the minimization problem is 

'■T ^ ^ f-T 



min E 



[xi + A5"xt + v^xi) dt 







z.X2coth(z.T) + ^tanhf — 
V V 2 



E 



E 



{{xi?+\s'^xi+v\x:f)dt 











(13) 



For 1/ = 0, the unique optimal strategy is given by (10) and the value of the minimization 
problem is 



min E 



[if + XS^xt) dt 



= E 









((i°)2 + A5?x?) dt 



(14) 



1 1 

— + -XSqT E 

T 2 16 



T 



{S?{T-t))'dt 



The strategies (9) and (10) are independent of the particular law of whenever 5° 
is a martingale. This has the immediate consequence that these strategies also minimize 
the following robust cost functionals. 

Corollary 1. Let Q be any set of equivalent probability measures Q on {Q,J-) under 
which the stochastic process is a square-integrable martingale satisfying (11) and for 
which {^l,J^, {J-t),Q) satisfies the usual conditions. When v = A7/?7 > 0, the strategy (9) 
minimizes the cost functional 



sup Eq 

and the minimal cost is given by 



[xt + XSiXt + v^xf) dt 



(15) 





i/X^ coth(i/T) + tanh f — - i inf Er 



min supEo / (±1 + XS?xt + v'^x'l) dt = sup Eq / {ixtf + XSUt + v"^ ix*,f) dt 
xex{T,x)Q<zQ L 7o J QeQ ^ Jo 

rT . 

, , r, nn ^ ' ' ' 

\ 2 J Au^ QgS 



Jo 



g^anh( ''^^, ^^ ^^dt 



An analogous statement holds in case u = 0. 



Remark 1. Let us summarize the positive and negative properties of the optimal strate- 
gies (9) and (10). 
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+ As shown by Theorem 1, the optimal strategies have a remarkable robustness prop- 
erty: they minimize the cost functional (7) (or its robust version (15)) for every 
martingale S^, regardless of the specific law of that process. If one accepts the cost 
criterion (7) and the assumption that the unaffected price process is well-described 
by martingale dynamics, then the strategies (9) and (10) will be optimal even in a 
situation of model uncertainty. These strategies are therefore robust with respect to 
model risk. 

+ The strategies (9) and (10) are given in explicit form and can be very easily im- 
plemented in practice, without the need for complex numerical methods. Their 
dependence on asset prices and model parameters is completely transparent. Also 
in this sense these strategies are very robust. 

- As a disadvantage, it should be noted that the strategies (9) and (10) can become 
negative. In practice, this will just lead to the early termination of the strategy and 
to the early liquidation of the asset position. But, as a consequence, the strategy 
will lose its optimality property in such a scenario. As discussed in (Gatheral & 
Schied 2011, Section 4), the probability that strategies become negative will be very 
small with reasonable parameter choices. The possible negativity of strategies can 
thus be seen as an effect that may be as negligible as the possible negativity of the 
unaffected price process in the Bachelier model, which is frequently employed in 
order execution. 



Theorem 1 is in fact a corollary of our next result, which applies to the situation of a 
general semimartingale S^. In this case, the optimal strategy can be conveniently written 
in terms of stochastic integrals with respect to the semimartingale 



Yt := --{Sf - So) + X f SUs. 
V Jo 

To be precise, stochastic integrals starting in t > will be defined as 

r dYu := r e« dYu = C dYu - f in dY^, 
Jt Jt+ Jo Jo 



and stochastic processes of conditional expectations such as 

l-T 



E 



sinh(z^(T — u)) dYu 



< t < T, 



(16) 



will be understood as the right-continuous version of this process, which exists since our 
underlying probability space satisfies the usual conditions. 



Theorem 2. For a general semimartingale satisfying (11), there exists a unique strat- 
egy in X{T,X) that minimizes the functional E[C(x) + Jq Xxt{S^ + jxt) dt] within the 
class X{T,X). For v = Xj/r] > 0, this strategy is given by 

x*t = sinh {u{T - t)) 



X 1 E [ j; sinh(i/(r - n)) dYu \ Fs 



sinh (uF) 



smh(i^(F - s))y 



ds 



(17) 
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and the value of the minimization problem is 



min E 

xeX{T,X) 

1 v2 



T 



C{x)+ / Xxt{Sl' + jxt)dt 





E 



X 



^X-'-XSo + v{^X^cothiuT)+ . 
2 V smn(i/i ) 



C{x*)+ I Xx^S^ + ^x^dt 


T 

E / smh{i^{T -u))dYu (18) 



T 



E[ sinh(z.(r-u))dy„|7-t]x2 



/g \ sinh(z^(T — t)) 

For = 0, the optimal strategy is 



dt 



x 







T-t 



T V 2./o (T-s)2 



E 



(r - u) dYu 



ds 



and the value of the optimization problem is given by 



dtriin E 

x(^XiT,X) 



C{x)+ / Axt(5° + 7Xt)(it 

JO 

^ . rT 



E 



T 



(19) 



-xSo + ??(— + -E / (r-^z)dy„ 



C(xO)+ / AxO(SO + 7x?) 

E[ /,^(r-^z)dy,|7-t]x2 







T-t 



In the proof of Theorem 2 it is actually shown that the respective strategies (17) or 
(19) minimize any cost functional of the form 



E 



T rT 

xt dYt + I {±1 + v'^xj) dt 
L Jo JO 



where Y is an arbitrary semimartingale with 

E[( sup \Yt\f] < oo. 

0<t<T 

The specific form (16) of y is not needed in the proof. Taking Y := ^S^ and u = thus 
yields the following corollary. 

Corollary 2. For a general semimartingale satisfying (11), there exists a unique 
strategy in X(T,X) that minimizes the expected costs E[C(x)] within the class X{T,X). 
This strategy is given by 



X 







T-t^^^ 1 



T 



rE 



T y 2r,J^ {T-sY 
and the value of the optimization problem is given by 



{T-u)dSl 



ds 



(20) 



min E 

x€X{T,X) 



C{x) 



E 



- XSo 



+r/— + -E 



T 



{T-u)dSl 



i-E 



m^{T-u)dSl\Ft]-.2 



T-t 



dt 
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Let = So + AI + A he the decomposition of the semimartingale into a local 
martingale M and an adapted process A of locally finite variation. Let us assume for 
simplicity that A has integrable total variation over [0,T]. Then (11) implies that M is 
uniformly integrable over [0, T] and thus a true martingale. Hence, the optimal strategy 
(20) can be written as 



T 



(T - u) dAu 



ds 



Here the drift A enters the optimal strategy basically in integrated form, and so one 
can expect that possible misspecifications of the drift may average out. This relatively 
stable behavior is in stark contrast to the direct dependence of optimal strategies on the 
derivative of the drift in models with transient price impact as found in Lorenz & Schied 
(2012). 

We conclude this section with the following outlook on possible generalizations of our 
problem. 

Remark 2 (Ambiguity with respect to market impact parameters). Our approach yields 
robustness with respect to the law of but it requires that the exact values of the market 
impact parameters, 7 and 77, are known. In reality, these parameters will be uncertain 
to some degree as well, but our results say nothing about the robustness with respect 
to these parameters. Also, our results rely in an essential way on the assumption that 
temporary impact is linear. It is not known to the author how the optimization problem 
from Gatheral & Schied (2011) can be solved for nonlinear temporary impact. Nonlinear 
price impact is often observed in data of financial transactions; see, e.g., Almgren et al. 
(2005). It also arises in a natural way when one attempts to solve our minimization 
problem under the additional 'no-buy' constraint X( < for t € [0, T]. 

Remark 3 (Multiplicative market impact model) . One of the shortcomings of the Almgren- 
Cliriss model is that asset prices can become negative when price impact gets too large. 
To avoid negative prices, Bertsimas & Lo (1998) proposed the following multiplicative 
price impact model, 

St ■= S^exp (jxt + r]xt), 

where 7, ?7 > 0. Naturally, this model goes well along with geometric Brownian motion 
as unaffected price process see also Forsyth (2011) and Forsyth et al. (2009). When 
defining C{x) = Jq Sfxt dt, one is led to the minimization of the functional 



E 







C{x) + A / Sfxt dt . (21) 



When 5" is a martingale, this cost functional can be expressed as 



E 



St / exp {-jxt + rjxt) {xt + Xxt) dt 
Jo 



Its minimization therefore boils down to the minimization of the classical 'action func- 
tional' Jq exp (7x4 + r]Xt) {xt + Xxt) dt over deterministic strategies x G (T, X); see also 
Remark 4.1 in Gatheral Sz Schied (2011). When S^ is not a martingale, however, the 
solution to the problem of minimizing (21) is not known to the author. ^ 
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4 Heuristic derivation of the solution 



The solution to our optimization problems is obtained by guessing the formulas (17) and 
(18) for the optimal strategy and value of the minimization problem and by applying 
stochastic verification arguments to show that these guesses are correct. In this section, 
we explain how the formulas (17) and (18) can be guessed heuristically. To this end, we 
assume that S" := S'' is a diffusion process with dynamics 

dSt = aiSt) dWt + h{St) dt, 

with sufficiently bounded and regular coefficients a{-) and b{-). Then, (2) implies that 

/ xt dSt + 11 Xf dt 
Jo Jo 

r-T r-T r-T 

--So- xta{St)dWt- / xtb{St)dt + i] / xfdt. 
^ Jo Jo Jo 

Hence, if a; € (T, X) is sufficiently bounded, 



C{x) 



^X^ - XSo 
2 

^-X^ - XSf) 



E 



C{x)+ I \xt[S^t + -fxt) dt =1x'^ - XSo + rfE 



■^i + xt{\St-h{St)) + u^xl)dt 



where b{x) = rj ^b{x). We thus define the value function 



C{T,X,S) := infE 



T 



f + xt{XSt-b{St)) + iy^xf]dt 



where the infimum is taken over all strategies x G X[T,X) for which Xt(7{St) dWt 
is a martingale (by (27) below this is actually the case for all x € X{X, T) as soon as 
cr(-) is bounded). Now we parameterize these strategies x € X{T,X) by their derivative, 
vt = Xt, and let V{T,X) denote the corresponding set of controls. For v € V{T,X) we 
define x^ := X + £ Vs ds. Then C{T, X, S) can be written as 



T 



+ x-t{XSt-b{St))+i^\x-,f) dt 



C(T, X, S) = inf E 
vev{T,x) 

Dynamic programming suggests that 

C{T - t,xl St) + (vf ds + xl{XSs -b{Ss)) + i^^«)2) ds 

should be a submartingale for every v G V(T, X) and a martingale as soon as v is optimal. 
Let us assume that C is smooth. Then an application of Ito's formula suggests that C 
should solve the following degenerate quasilinear PDE: 



Ct = la\S)Css + biS)Cs + X{XS - b{S)) + u^X^ + m{{v^ + vCx). 

rT 



(22) 



In addition, the fuel constraint vt dt = —X required from strategies in V(T, X) suggests 

(23) 



that the value function C should satisfy a singular initial condition of the form 

ifX = 



limC(r, X, S) 

no 



+00 iiX^O. 
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The intuitive explanation for this initial condition is that a nonzero asset position with 
no time left for its liquidation means that the liquidation constraint has been violated. 
Excluding this violation requires an infinite penalty; see also Schied et al. (2010) for similar 
effects in utility maximization for order execution. 
To solve (22), (23), we make the ansatz 

C(T, X, S) = uX^ coth(z/r) + Xh{T, S) + g{T, S). 



Then 



2 v2 



Ct 

Cx 
Cs 
Css 



v'X 



(sinh(z^r))2 
2vX coth(z^r) + /i(r, S) 
Xhs{T,S) + gs{T,S) 
Xhss{T,S) + gss{T,S). 



+ XhTiT,S)+gT{T,S) 



Plugging this into (22) yields 



(sinh(i/T))2 



+ XhT{T,S) + gT{T, S) 



1 



Xhss {T,S) + gss (T, 5) + bXhs {T,S) + bgs (T, S) 



+X{XS - b{S)) + 1/2^2 - ^ ( 2iyX coth(z^T) + h{T, S) 



Since 



we obtain 



■^X^ - ^ (2uX coth(i/r) + h{T, S) 

7 ■'ff' 2 - lyXcotHuTMT^s) - h{T,sf, 

(smh(i/i jj^ 4 



XhT{T,S)+gTiT,S) 

= ^a^ [xhssiT, S) + gss{T, 5)) + bXhs{T, S) + bgs{T, S) 

+X{XS -b{S)) - uX coth{uT)h{T, S) - -h{T, Sf. 



(24) 



Equating all terms containing X yields 

Ht = ]^a'^hss + bhs + XS -b{S) - ucoth{uT)h. 

As initial condition we take h{0, S) = 0. This initial-value problem can be solved by 
means of the Feynman-Kac formula (Karatzas & Shreve 1991, Theorem 5.7.6). To this 
end, we define h{t,x) := h{T — t,x) for some fixed T > 0. Then h satisfies the terminal 
condition h{T, x) = and the PDE 



-ht + ucoth{u{T - t))h = -a^hss + bhs + XS - b{S). 
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It follows from the quoted Feynman-Kac formula that, for t < T and Y is as in (16), 
rT 



h{t, x) 



[XSu - b{Su)) exp 



exp 



sinh(i/(r - t)) 



E 



V coth(i/(T — r)) dr ) dl^ 

T 



V coth(;/(T — r)) dr ) du 
St = x 



St = X 



sinh(z^(T — u)) dY^ 



St = X 



Hence, 



h{T -t,St) = h{t:St) 



1 



sinh(z^(r - t)) 



E 



and in particular 

h{T,So) = h{^,So 



1 



-E 



sinh(i/(T — u)) dYu I J-'t 



sinh(z^(T — u)) dY^ 



sinh(z^T) 

The PDE for the function g is obtained by taking X = in (24): 

9T = ^cr'^gss + bgs - ^/i(T, Sf. 
The initial condition must be ^(0, S) = 0. It follows that 

g{T,So) = -E 



^{h{T-t,St))^dt 



1. 



-E 



E[f^smh{i^iT-u))dY^\Tt 



sinh(i/(T-t)) 



dt 



Putting everything together, we obtain 



c(r,x,5o 



i^X"^ coth(i/r) + 



X 



1e 

4 



sinh(z^T) 



E 



sinh(z^(r-n))(iFn 



^ ^E[ sinh(zy(r - u)) dYu\J^t^^^ 



dt 



(25) 



sinh(i/(T - t)) 

Note that this formula is independent of the particular characteristics a{-) and b{-) of 
the dynamics of S. We therefore may guess that the formula remains true for even more 
general semimartingales, which lead us to asserting (18). 

Standard arguments in control suggest that the optimal strategy x* is defined through 
that V* that attains the infimum in (22). More precisely, x* should be the solution of the 
ODE 



1 



Cx(T - t,xl,St) 



-vxl coth(z^(r - t)) - hi{T - t, St). 



X 



This ODE is solved by 

x'l = sinh ((T — t)z^) 

— sinh ((T — t)v) , , , 

' ^ [sinh (i/T) 2 Jo 

and this formula suggested the assertion (17). 



1 /■* h{T-s,Ss 



■ ds 



sinh (i^T) 2 Jq sinh [i'{T — s)) 

X 1 /■* E [ /J sinh(z/(r - u)) dYu \ J^s 



( sinh {i^{T — s)))' 



ds 
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5 Proofs 



Proof of Proposition 1: Let x € X{T,X) be given. We have 

r r r 

C{x) = I S^xt dt = S^xt dt + r] xfdt + j {xt — xo)xt dt. 
Jo Jo Jo Jo 

The rightmost integral is equal to 



JXq + 1 ^txt dt = ^xl = ^X"^. 



Hence, 



C{x)+ I XxtiS^ + ^xt) dt 
2 



^2 + ^ Sf±tdt + r](^j^ xfdt + j^ {XS?xt + iy^xf)dty 



(26) 



The first two integrals on the right-hand side of (26) belong to L^{F) by (11) and (12). 



Next, the fuel constraint xt dt = —X and Jensen's inequality imply that 



Xt 



J Xs dsj <{T-t) J xldse 



(27) 



It follows that also the rightmost integral in (26), and in turn the entire expression (26), 
belong to L^{F). 

Now suppose that S'^ is a martingale. When J xt dS^ is a true martingale, our result 
follows by applying integration by parts. The martingale property of J xtdS^, however, 
is not clear without additional integrability conditions. We therefore proceed as follows. 
Since we have assumed that the filtered probability space (17, F, (J-j), P) satisfies the usual 
conditions, the rightcontinuous martingale (5^) is optional; see Theorem 65 in Chapter 
IV of Dellacherie & Meyer (1978). Thus, the process {S^f^j^) is equal to the optional 
projection of the constant process t ^ Sj.. Hence, 





r /"^ 1 




rp 




rp 


E 


/ S^±t dt 
L Jo J 


= E 


[ SUxt 


= E 


' / S^rdxt 

L Jo J 



S^dxt =nST{xT-xo)] = -SoX, 



where in the second step we have used Theorem 57 in Chapter VI of Dellacherie & Meyer 
(1982) and the fact that \SiXt \ dt belongs to L^(P) due to our assumptions on and 



x. This concludes the proof. 



□ 



Proof of Theorem 2: Let x E <Y(T, X) be given. Integrating by parts yields 

/ SgXsds + r] / XSgXgds = S^xt- SqXq + t] / XgdYs- 
Jo Jo Jo 

Hence we get from (26) that 

7^2 



(28) 



C{x) + [ XxtiS? + 
Jo 



7Xt) dt = jX' -SoX + r,i^ 1^ xt dYt + {it + v'xt) dt ) . 
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Our problem therefore reduces to minimizing the expression 



E 



T 



T 



xtdYt+ / {xi + v^xi)dt 



(29) 



over X S X{T, X). We first consider the case > and show that the optimal strategy is 
(17) and that the minimal value of the cost functional (29) is 



vX"^ coth(z^r) + 



X 



smh{vT) 
T TT^r rT 



E 



E[ // sinh(i/(r - u)) dYu\Ft\\2 



sinh(z^(r - t)) 



sinh(z^(r — u)) dYu 
dt 



(30) 



Let us study the various expressions in (30). Integrating by parts and using 10 = 0, 
we find that 

I smh{u {T - u)) dYu = sinh{u {T - t))Yt + u [ cosh{u {T - u))Yudu 
Jo Jo 



and hence that 
rT 



[ smli{ij{T -u))dYu = -smh{ij{T -t))Yt + iJ [ cosh{ij{T - u))Yudu. (31) 
Jt Jt 



The right-hand side is square-integrable due to (11). 
Let M be a rightcontinuous version of the martingale 



E 



sinh(z^(T — u)) dYu J~t 



i/E 



cosh (z/(T — u))Yu du 



0<t<T, 



which exists since the underlying probability space was assumed to satisfy the usual 
conditions. Clearly, M is a square-integrable martingale. 
Then 



1 



sinh(z^(r - t)) 
1 



E 



and so 



sinh(z^(r - t)) 
dHt = i/coth(i/(r - t))Ht dt + 



sinh(z>(T — n)) dYu \ Tt 

Mt- [ sinhiu {T - u)) dYu 
Jo 



1 



sinh(z^(r - t)) 



dMf - dYf. 



(32) 



Moreover, by (31), 

Ht = -Yt + 

and so 



sinh(z^(r - t)) 



cosh(z^(T — n))E[y„ | T ] du, 



\Ht\ < \Yt\+E[ sup \Yu\\Tt]. 

t<u<T 



We thus conclude from (11) that 



sup \Ht\ G L 
0<t<T 



(33) 



(34) 
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Therefore, we may define as a rightcontinuous version of tlie martingale 

1. 



r 




/ H^du 




. Jo 





0<t<T. 



We also define 



Hldu 



Nt + - I Hidu. 



It follows that 



dGt = dNt + -,Hf dt. 



1 

4' 



(35) 



Now let X € X(T, X) be given and define 

Ct := i^xf coth(z^(r - t)) + xtHt + Gf 

When the theorem is correct, Gt should be the minimal cost (29) for liquidating the 
position xt over [t,T]. Hence 

ft ft 

^2 , ,.2 2\ 



Ct:= / XsdYs+ / {±i + u''xi)ds + Gt 
Jo Jo 

should be a submartingale for any x S X(T,X) and a martingale for the optimal x*. 
Let us first prove the following claim: 



sup \Gt\eL\¥) 

0<t<T 

First, it follows from (28), (11), and (27) that 



for ah xeX{T,X). 



(36) 



sup 

0<t<T 

From (12) and (27) we get next that 



(37) 



{x'i + i^x1)ds G L^(P). 



Furthermore, the fact that the function t i— )■ t cotht is bounded on [0, z^T] by a constant c 
implies together with (27) that 



0<t<T 

Using again (27) and now (34) gives us 



sup uxf coth(z^(r -t)) <c I ±1 dt € L^(P). 
<t<T Jo 



Finally, 



sup \xtHt\ G L\¥). 

0<t<T 



sup \Gt\e L\¥) 

0<t<T 
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follows from (34) and Doob's L^-maximal inequality. Our claim (36) now follows by 
putting everything together. 

Using (32) and (35) we compute 



dCt 



Xt dYt + 



2^2 



Xf- + v Xf- + 2xtXti' coth{i'{T — t)) + 



I' X 



-xti/ coth(i/(r - t))Ht + -Hi 



dt — Xt dYt + 



(sinh(i/(T - t)))2 

Xt 



sinh(z>(r - t)) 



+ ±tHt 
dMt + dNt 



Xt + xtu coth(i/(r - t)) + -Ht dt + 



Xt 



sinh(z/(r - t)) 



dMt + dNt. 



(38) 



Hence, C is the sum of a nondecreasing process and a local martingale. 
Let us introduce the stopping times 



Tn ■■= inf jo < t < T \xt\ > n| A T, 



with the convention inf = +oo. Then 







sinh(i/(r — s)) 



dM, 



is a true martingale for < t < T and n € N. Therefore E[CfAr„ ] > Co = Cq. By (36) 
we may pass to the limit t f T and n f oo and obtain 



E 



xt dYt + / {xl + u'^x'i) dt] = E [ C7t ] >Cq. 



(39) 



Note that Co is equal to the asserted minimal value (30) of our optimization problem. 

Since E[CtAr„ ] is nondecreasing in t and n, we have an equality in (39) if and only if 
E[CtAr„ ] = Co for all t € [0, T) and n G N, which by (38) holds if and only if x satisfies 
the ODE 



Xt = -xtv coih.{v{T - t)) - ^Ht, 



■^-a.s. for almost every t £ (0, T). (40) 



When a solution to this equation exists in X{T,X), it will be unique in this class by 
standard arguments. Thus, there is at most one optimal strategy. 
Note that the strategy (17) can be written as 



sinh {iy{T - t)) 



X 



1 



sinh {fT^ 2 Jq sinh (z^(T — s)) 



ds 



and hence solves (40). We show now that x* € X{T,X). First, we clearly have Xq = X. 
We show next that x^ as t '[ T. We have 



1 



/ 7 — : ds = — log 

Jo sinh (i/(T - s)) u 



1 /sinh {u{T - t)) cosh(i/r) - 1 

V sinh(i/r) cosh (i/(T - t)) - 1 ) ' 



As t ^ T, the right-hand side behaves like 



1 , cosh(z^r) -11, 1 

- log ■ + - log — — -r, 

V sinh(i/i j V lv\l — t) 



(41) 



(42) 
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This function belongs to L'^[0,T], and with (34) we get -i' as t ^ T. We finahy 
have to show that dt G L^(P). We know that x* satisfies the ODE (40) and that 

Jq Hfdt e Li(P). Next, 



\xl coth(i/(T — t))| < ci + C2 / 

Jo 



ds ■ sup \Hs 



sinh [u{T - s)) o<s<T 
for certain constants ci,C2 > 0. From (41), (42), and the fact that 



log 



1 \k 



l-t 



dt = k\ 



for A; G N, we thus conclude that 

/ ix*tf dt < C4 + C3 • sup \Hs\^ G L^(P). 

Jo 0<s<T 

Hence x* G X{T,X), and the proof for the case 7 > is complete. 

The case 7 = can be analyzed by passing to the limit 4, 0. To this end, we note 
first that for < t < u < T 

sinh(z^(T — n)) , , sinh(z/(T — n)) T — u 

0< -V-^— ff < 1 and lim- ^ ^ 



sinh(i/(T - t)) 



L/iO sinh(z^(r -t)) T -t' 



It hence follows from the dominated convergence theorem for stochastic integrals (Protter 
2004, p. 267) and (Protter 2004, Chapter V, Theorem 2) that 



sup 

0<t<T 



sinh(i/(r - u)) dY^ J^{T - u) dYu 



sinh(z>(r - t)) 



T-t 



in L^(P) as iO. 



Therefore, letting 

Xi := sinh (i^(r — t)) 



X 1 /■* E [ sinh(z^(r - u)) dYu \ J^s 



sinh (lyT) 



sinh {h'{T — s)))' 



ds 



we have supQ^^^j^ \x^ — | — t- in L^(P). Passing to the limit J, in (40) leads to the 
ODE 



xt 



-xt 



T-t 2{T -t) 



E 



(T - u) dYu 



P-a.s. for almost every t G (0, T), 

(43) 

which is solved by x^ . Hence, we also have x\ x\ in dt ® dP- measure. It follows that 



lim inf E 

v\a L Jo 



T 



T 



\dY,^ \ [[±\f^v\x^,f)dt 



> E 



T 



xUYt+ I {x^tfdt . (44) 



For each fixed > 0, the expectation on the left-hand side is given by (30), and the latter 
expressions converge to 



X 

— + — E 
T T 



(T - u) dYu 











4 





E[ //(T-n) dYu\Ft 
T-t 



dt 



(45) 
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which is hence larger than or equal to E[ dYt + /g^(x^)^ dt] . Finally, the optimality 

of follows from (44) and the facts that each x'^ is optimal and that 



E 



xt dYt + 



fx? + i/'^Xf) dt 



— > E 









Xt dYt + 



^0 



dt 



as V 10 for each x e X{T, X). 

Proof of Theorem 1: When is a martingale and > 0, 



E 



T 



sinh(i/(T — u)) dY^ 



'- Jt 



Tt 



AE 



sm\y{v{T - u))Sldu 



Tt 



XSf smh{i^{T -u))du 
XSt- sinh(z^(r - t)) tanh 



□ 



Thus, (30) reduces to (13) in the martingale case. Moreover, (17) reduces to (9), because 



tanh (^) 



1 



sinh(z/(T - t)) 1 + cosh(i/(T - t)) ' 
This concludes the proof of Theorem 1. 



□ 
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